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I. INTRODUCTION 



W^'ip^q) = J d'xe^'^-^p,S\ [J^ix),rm\p,S), (1) 
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We discuss the determination of deep-inelastic hadron structure in lattice QCD. By using a 
fictitious heavy quark, direct calculations of the Compton scattering tensor can be performed in 
Euclidean space that allow the extraction of the moments of structure functions. This overcomes 
issues of operator mixing and renormalisation that have so far prohibited lattice computations of 
higher moments. This approach is especially suitable for the study of the twist-two contributions to 
isovector quark distributions, which is practical with current computing resources. While we focus on 
the isovector unpolarised distribution, our method is equally applicable to other quark distributions 
J^^", and to generalised parton distributions. By looking at matrix elements such as {7r*|T[y^(x)^''(0)]|0) 

• (where and A" are vector and axial-vector heavy-light currents) within the same formalism, 

' moments of meson distribution amplitudes can also be extracted. 

(N ' 
o 

Q, 

Lattice QCD offers the prospect of exploring the structure functions probed in deeply inelastic scattering (DIS) and 
' other high-energy experiments from first principles. By comparing to accurate experimental data, such calculations 
provide stringent tests of QCD. They also allow the extraction of information on hadron structure which is not 
. currently available from experiment, e.g., the transversity distribution 5q{x). The structure functions describe the 
^ ' hadronic part of the DIS process, viz., the hadronic tensor 
l> ■ 

o 
o 

. 

where p and S are the momentum and spin of the external state, q is the momentum transfer between the lepton and 
■ the hadron, and is the electromagnetic current. Using the optical theorem, Wg'^ can be related to the imaginary 
\ part of the forward Compton scattering tensor 
> . 

^ : = / d^xe'i-^p,s\T[j^{x)rm \p,s). (2) 

Since lattice QCD is necessarily formulated in Euclidean space, direct calculation of the structure functions is 
' challenging because of the analytical continuation to Minkowski space that is required 0, Q- In addition, such 
^ ,. a calculation would involve all-to-all light-quark propagators, and is therefore numerically demanding. For these 
\^ ■ reasons, beginning with pioneering works of Refs. lattice studies of the deep-inelastic structure of hadrons have 

I ^S focused on calculations of matrix elements of local operators that arise from the light-cone operator product expansion 
(OPE) d H 0, i, of the currents 

_ _ I 

T[J^{x)rm^Y.^,{x\i?) x,,...x,^Or'^-^-{^), (3) 

i.n 

where the Ci are the perturbatively calculable Wilson cocfhcicnts that incorporate the short-distance physics, and the 
sum is over all local operators, ^^'^^i - '^" with the correct symmetries. ^ is the renormalisation scale. This expansion 
enables the investigation of Tg'^ via the knowledge of hadronic matrix elements of local operators. The analytical 
continuation of these matrix elements from Euclidean space to Minkowski space is straightforward. However, a number 
of difficulties arise in this approach because of the lattice regularisation. Firstly, the non-zero lattice spacing breaks the 
symmetry group of Euclidean space-time from 0(4) to the discrete hyper-cubic subgroup H(A), consequently modifying 
the transformation properties of the local operators in the OPE. In general, operators belonging to different irreducible 
representations of 0(4), which span the right-hand side of the OPE in Eq. lO, mix unavoidably in the lattice theory 
since i?(4) has only a finite set of irreducible representations. For twist- two (twist = dimension - spin) contributions, 
this becomes particularly severe for operators of spin n > 4 as they mix with lower dimensional operators and the 
mixing coefficients contain power divergences. Currently this restricts the available lattice calculations t o op erators 
of spin n — 1,2,3,4. For higher-twist operators, such power divergences are generally unavoidable llOl, [ill [lj|. A 
second issue is that the matching of the lattice regularisation to continuum renormalisation schemes [l3l Il4l Il5l , 
in which the Wilson coefhcients are calculated, becomes more involved as n increases. 
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In this paper, we discuss an approach to determining matrix elements of higher-spin, twist-two operators in Eq. 
This approach is based upon directly studying the OPE on the lattice, as was first investigated in kaon physics in 
Ref. 17]. A similar technique has also been applied to determine Wilson coefficients non-perturbatively jl^l and 
extract the lowest moment of the isovector twist-two quark distribution '19*1 (our method is related to this latter 
work but improves on it in a number of ways). In our proposal, one simulates the Compton scattering tensor 
using lattice QCD, with currents coupling the physical light quarks, ip{x), present in the hadron to a non-dynamical 
(purely valence), unphysically heavy quark, ^(x).^ The introduction of this heavy quark significantly simplifies 
the calculation of isovector matrix elements because it removes the requirement of all-to-all propagators. After 
performing an extrapolation to the continuum limit, the lattice data for the Compton tensor are compared to the 
predictions of the OPE in Euclidean space to extract the matrix elements of local operators in Eq. Q, directly in 
the continuum renormalisation scheme in which the Wilson coefhcients are calculated. This approach also removes 
the power divergences, thereby enabling extraction of matrix elements of higher spin (n > 4) operators for twist-two 
operators with a simple renormalisation procedure. These matrix elements determine the Mellin moments of the 
structure functions which are identical in Euclidean space and Minkowski space and their analytical continuation is 
trivial. Finally, the chiral and infinite volume extrapolations can now be performed at the level of the local matrix 
elements using chiral perturbation theory j21j, ,2^ 2^, 2As> 2hs 2^ 2L 2^ JSJ • 

The matrix elements obtained via the above procedure are completely independent of the mass of the unphysical, 
heavy quark and are indeed physical quantities. This is because such a quark can only propagate between the 
bilocal currents, and the OPE relegates its short-distance information to the Wilson coefficients. In addition to the 
numerical advantage, it also proves useful to introduce a fictitious heavy quark for other reasons. Firstly, the presence 
of the heavy scale suppresses long distance correlations between the currents in a similar way to a large Euclidean 
momentum. Combining both the heavy quark mass, m^, and momentum injection, q, at the current allows us to 
control the behaviour of the OPE precisely at moderate and q^. The only constraint is 



where a is the coarsest lattice spacing used in the calculation. Secondly, the non-dynamical nature of the heavy quark 
automatically removes many contributions (for example, so-called "cat's ears" diagrams - see Fig. PTd) below) that 
are higher- twist contaminations in traditional DIS. 

In Section^ we review the formalism of DIS with heavy quarks before discussing the extraction of the moments of 
twist- two parton distributions from lattice correlators in Section IIIII Finally in Section IIVI we broaden the analysis 
to investigate moments of meson distribution amplitudes. 

II. FLAVOUR CHANGING CURRENTS AND HEAVY QUARKS IN LEPTON-HADRON 

DEEP-INELASTIC SCATTERING 

The roles of quark and hadron masses in deep-inelastic scattering have been well studied. Target mass effects were 
first discussed by Nachtmann [2OI and extensively investigated throughout the 1970s, following the observation of 
the precocious scaling of the structure functions |^ [s^]. Away from the Bjorken limit, they result in significant 
contributions which arise from the OPE being an expansion in terms of operators belonging to definite irreducible 
representations of the Lorentz group. These contributions scale as powers of M^/Q^, where M is the target mass and 
= —q^, and can be summed exact lv jSOl l33l l34L l35 | . The effects of the struck and produced quark masses were 
also comprehensively investigated [s^, l34l3()|. These target and quark mass effects lead to ^ scaling [s^, 113, 0, HI] , 
and are particularly relevant at moderate values of Q^. Since currently available lattice cut-offs are 1/a ~ 3 GeV, 
it is important to include these mass effects in the application of the OPE on the lattice, because of the condition 
in Eq. In this section we present the OPE in Euclidean space relevant for computing higher moments of parton 
distributions on the lattice with these mass effects taken into account. 

We consider fictitious currents that couple light up and down quarks to unphysical heavy quarks of mass rriif . We 
focus on a purely vector coupling, leaving the discussion of other possible currents to the end of the section. We define 




(4) 



J^,^^) = *(x)7^^(a;) + ^(x)7^*(x) , 



(5) 



Such fictitious currents have been used to study quark-hadron duahty in heavy quark effective theory I20I . However, in this context the 
heavy quark is not a valence quark. 





FIG. 1: Contributions to the Compton scattering tensor. Diagrams (a), (b) and (c) correspond to the leading twist contributions. 
Diagram (c) (the "box diagram") involves gluonic operators and vanishes for the isovector combination, Eq. JTJ. Diagram (d) 
(the "cat's ears diagram") is higher-twist and absent in our analysis. Diagram (e) includes leading- and higher-twist terms and 
is discussed in the main text. The thick lines correspond to the heavy-quark propagators, the shaded circles to the heavy-light 
currents and the large shaded regions to the various parton distributions. 



and construct the Euclidean Compton scattering tensor 



E 



(6) 



oo, T^"^ is given by the leading-twist contribution, the "handbag diagrams" in 



(henceforth all momenta arc Euclidean) 
In the limit ^ oo or 

Figs, m (a) and (b). The "box diagram", Fig. (c)^, which involves purely gluonic operators after the OPE, is 
strongly suppressed in our approach and is completely absent in the study of the OPE of the isovector Compton 
scattering tensor 



rpflV 



(7) 



This makes the extraction of moments of the isovector quark distributions practical, and we focus on this case in this 
paper. 

At moderate and m>f , higher-twist terms also contribute. However, the non-dynamical nature of the fictitious 
heavy quark entirely eliminates the higher-twist contributions involving more than one quark propagator between the 
currents, e.g., the "cat's ears diagram" in Fig.n(d). The diagrams in Fig.n(e) contain pieces that contribute to the 
twist- two operators in Eqs. \i'2.\ and H13|l. and also higher- twist terms that are discussed below. 

The twist-two contributions to the OPE in T^''^ are from 



[iD + g)2 + m% 



(8) 



Dfj- _ JJ^^ 1 are included 



and a similar term. Fig. ^(b), in which fj. and q —q. The derivatives [D^^ — i 

to account for the soft transverse momentum of the struck quark; they are covariant in order to maintain gauge- 
invariance. 



^ In Fig. 0(c), we specify that the large momentum, q^, flows through the three light-quark lines; the contributions in which these quarks 
have soft momenta are already included in Figs.Q(a) and (b). In principle, these gluonic contributions can be disentangled from their 
different behaviour. 
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The OPE of the bilocal currents in Eq. © is now given by a Taylor expansion of the propagator in the above 
expression, 



771 xj. 



771^! 



-2i q-D 



[iD + qY + 77i\ 



£,2 



n=0 



(9) 



where = —q^. The pole mass of the heavy quark is not a measurable quantity and we replace it with the mass 
of a meson composed of the unphysical heavy quark and a light quark, Mip = m^, + ^a, where a is the binding 
energy and is of ©(Aqcd)- This meson mass can be directly computed on the lattice. The parameter a can be 
extracted from experiment or lattice calculations [s^ . however we view it as an unknown to be determined in the 
procedure of studying the OPE on the lattice. The term in the denominators in Eq. ^ gives rise to higher-twist 
contributions, such as those arising from Fig.n(c), if we instead Taylor expand with respect to ^ ~q2 ^ ■ These 

higher-twist contributions scale as powers of ^}^QXili^^ + 1^%) ^^nd can be neglected for q^ or 7n\ sufficiently large. 
For more moderate scales, they may become important (particularly for large n) and can in principle be studied in 
our approach. However, this is beyond the scope of this work. Instead, we replace the denominator in Eq. Q by 



Ml 



aM<L 



(10) 



where the unknown parameter (3 represents terms of 0(Aqqj-,), including these higher- twist effects and sub-leading 
heavy-quark mass effects. In doing so we have neglected the n dependence of (5 since it is suppressed by powers of 
the strong coupling. To remove these higher- twist and mass uncertainties, one might consider using a fictitious heavy 
particle that does not interact strongly. However, issues such as the gauge dependence of the resulting Compton 
tensor would need to be investigated thoroughly and are not discussed here. 



After making this expansion and considering only the symmetric combination (m is the light quark mass) t^^^ 



combining the two leading contributions then leads to 



71+2- 



(29mi)---(29mJ 



even 



Q 



2n 



^qtJ,UfJ,i...ll„ 



g2^M^ ^ C: 



, (2g^J...(2g^J 



n=2 
even 



12n 



+2iS>"'{77i^ -to)^C„ 



Tl = 

even 



(2g^J^..(2g^J 

Q2n 



^4, 



(2g^J...(2g^^) ^^}^ 



odd 



V'7 



{M1 



and 



^r^^ll...^l■n 



jD'""M -0 - traces , 



traces , 



where we have introduced the operators (the braces indicate symmetrisation of the enclosed indices) 



(11) 



(12) 



(13) 



and the various perturbatively calculable Wilson coefficients C„, C,^, C^' and C„ depend on q^, /^^ (the renormalisation 
scale) and the heavy quark mass (38i. .391. liollZH . The first set of operators are the usual twist-two vector operators 
that enter into textbook DIS analyses. The second set are chiral-odd, twist-three operators whose matrix elements 
correspond to the parton distribution (x) [i^ l43l| . 

Taking spin-averaged hadron matrix elements of this expression then leads to 



=2 



nCW(r?)~27^Cr_^,(yy) 
(p- 7i(n — 1) 



(2) 



(p • qf 



77(71 — 1) 



(14) 



+ 4: 



p-q 



g2 i77-i)7jcl^lM-^v'c':Mv) ^„VM2) f , 

'92 n(n-l) ^n^^n^AV) 



o2/^(3) 



Q2 n{77 ~ 2)CW(77) ~ 277(2n - 3)c1^1M + Sy^C^J^Uv) / _ ^^Z^^') (^)) 

q^ 71(77 — 1) n y n ^ n-1 J 



even 
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where we have defined 



(M is the proton mass, = — Af^) and the hadronic matrix elements of the local operators in Eqs. (I12|l and H13() as 

Y,{v.S\0^^'-^-\v.S) = -traces] , (16) 

s 

Y,l^,S\d%'-^-\v.S) = iMA5(A.2)[pM,^^pM„ .traces] , (17) 
s 

(the and are dimensionless and real). In Eq. p4|l the Cn^\r]) are Gegenbauer polynomials that arise from 
the trace subtractions in Eqs. lfTB|l and ((TTIl and account for the target mass effects [s^, 113, 13 H^- If we now 
choose the rest frame of the proton, p — (0, 0, 0, i M) and parameterise q = (0, 0, \J%—Q^, i qo), then the symmetric 
combination of {fi, — {3, 4} is 

oo 



n=2 
even 



where 



^"Q2 n(n-l) +^n^L.„_iW 



(19) 



Equation (|18|l is the central object studied in the following section where we will see that lattice calculations of it 
would allow extraction of the even moments ^^(/i^), since the Wilson coefficients and kinematic factors that determine 
f{n) are known. However we discuss some more general aspects of T^"^ first. 

Once we have determined the A'^{fj.'^), we can also use the diagonal elements of T^J^^ to extract information on the 
moments of the twist-three quark distributions e^{x) from the same vector- vector current correlator. Experimentally, 
this distri bution is difficult to extract, being measurable only in charged current DIS, semi- inclusive DIS or Drell-Yan 
processes 42, 43]. The only currently available determination uses data from the CLAS collaboration and requires 
assumptions about the fragmentation functions that enter semi-inclusive DIS |l45l] . Any information on the moments 
of 6^(2:) from lattice QCD (whether in our approach or by direct calculation) would be useful. Additionally, the odd-n 
moments ^^(/i^) and ^^(/i^) (which determine the valence combinations of quark distributions) can be extracted 
from a correlator of vector and axial- vector currents just as for ^3(0;) in neutrino scattering. 

Similar methods can be used to calculate moments of the helicity and transversity parton distributions. The 
moments of both the twist-two and twist-three helicity distributions and the twist-two transversity distribution may 
be determined from suitable antisymmetric pieces of T^''^. Since we are not restricted by physical scattering processes, 
we can also consider correlators of unphysical currents: e.g. scalar, pseudo-scalar and tensor. The analysis for such 
currents is again very similar to the case that we have discussed and by using such currents we can investigate 
distributions that are not experimentally accessible in DIS; for example the transversity distribution can be accessed 
independent of the quark masses by looking at the scalar-vector correlator. Since these currents are not conserved, 
they are no longer scale independent and the appropriate anomalous dimensions must be used [4^ |47| . Finally, 
the off-forward matrix elements of the various twist-two and twist-three operators that determine the moments of 
generalised parton distributions can also be studied through analysis of the off-forward Compton scattering amplitude. 



III. EXTRACTION OF MOMENTS FROM LATTICE CALCULATIONS 

With numerical investigation of T^^^^ {p,q) for varying Q^, qo and M^j,, Equation H18|) provides a means to extract 
the moments A'^(fi) for n > 4 without requiring power subtractions and complicated renormalisations. In order to 



(a) (b) (c) 

FIG. 2: Quark contractions in the four point correlator. Light shaded circles are the proton source and sink whilst the dark 
shaded circles are the heavy-light currents. The thick line indicates the heavy quark propagator. Diagram (c) is quark-line 
disconnected (as indicated by the representative gluons). 



calculate T^'^^^ (p, q) we consider the following four-point Euclidean correlator HIH (with r > 0) 

G|';)(p,q,t,r;r) = J] ^ e'P-e''i-yr^„(0|x„(x, i))7^,^(y -f z, r -f |)j;,^(z4)x/3(0, 0)|0) (20) 

x,z y 

^ ^ ^ e«(p-PN)-Xgj(pN-PN')-Zgjq-yg-(SN + Siv')|p^^ 

x,z y N,N' s.s' 

x(0|Xa(0)|£;Ar,PN,s)(£'Ar,PN,sk^^,,^(y,T)75^,^(0)|£;Ar',PNSs')(£:Ar,,pNSs'|x^(0)|0} 

^.Eot g.q.y ^ r0„(O|xa(O)|i?o, P, s){Eo, p, s|7^,^(y, r)j;,^(0)|i?o, P, s') {Eo,p, s'|X/5(0)|0) , 
y s.s' 

where x is a dimensionless interpolating field for the proton and ^ is the lattice version of Eq. (O . Here F is a 
Dirac matrix which can be chosen as r4 = i(l -I- 74) for the components of the amplitude we are considering and 
Eo{p) is the ground state energy. The sums on N and N' are over all possible eigenstates of the Hamiltonian. Then, 
defining the two-point Euclidean correlator as 

G(2)(p,t;r) =^e'P-^r^<,(0|xa(0)x/3(x,i)|0), (21) 

X 

we can determine the Compton amplitude from the Fourier transform of the ratio of these correlators: 



Gl^r\p,q,t,T;Ti) 

lim — — 

t^oo G(2)(p,i;r4) 



(22) 



where a is the lattice spacing. Since G^^^^(p, q, t, r; F) falls off with r over a characteristic time (M^a) ^, the Fourier 
transform in Eq. H22|l will be well approximated in practice. 

The above equations also hold if the heavy-light current is replaced by the usual light-light current. However, using 
a quenched heavy quark greatly simplifies the numerical work in extracting T^'^^ {p, q) without altering the non- 
perturbative physics, the Mellin moments, we are interested in. After performing the Wick contractions of the quark 
fields, G^^'^^ (p, q, i, t; F) is given by the three different arrangements of quark propagators shown in Fig. El where the 
thick and thin lines represent the heavy and light (physical) quarks respectively. For the isovector channel obtained 
via Eq. (|7|), diagram[21 (c) does not contribute. In the remaining two diagrams, the Wick contractions can be computed 
with the technique of extended propagators. Many values of the heavy quark mass can be studied with only a very 
small increase in computational cost. For a light-light current, we would need the additional Wick contractions shown 
in Fig. 13 and would require light all-to-all propagators even in the isovector channel. 

In order to avoid the problems of operator mixing discussed in the introduction, a continuum extrapolation of the 
Compton amplitude must be performed before it can be used to extract the moments A!^ (p) in Eq. H18|l . This requires 
calculations at a number of different lattice spacings. On the other hand, the chiral and infinite volume extrapolations 
can only be performed after we have extracted the local matrix elements since the current-current correlator is not a 
low-energy observable and cannot be treated reliably in effective field theory. Also, to use Eq. (|18|l . the heavy-light 
meson and proton masses must be extracted from appropriate two-point correlators. 

Assuming these tasks have been performed, we may now consider what is required to obtain useful information from 
these calculations. Ideally we would extract all integer moments from the lattice and be able to uniquely determine 





(a) 



(b) 




FIG. 3: Additional quark contractions present in the four point correlator for light-light currents but absent in the case we 
consider. These diagrams correspond to higher-twist contributions (diagram (d) in Fig.0. Light shaded circles are the proton 
source and sink whilst the dark shaded circles are the light-light currents. 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 



10 
n 



15 



20 



0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 







* * * 



10 
n 



15 



20 



FIG 



4: Behaviour of the function f{n) for two different parameter sets. On the left, 

l2 o or ry„\r2 



3.54 GeV, = 1.5 GeV^ and 

qo — 2.76 GeV while on the right we have chosen M* = 2.1 GeV, = —3.85 GeV^ and qo = 1.98 GeV. The stars and boxes 
correspond to q = 0.4, 1.2 GeV respectively and we have chosen M — 1.2 GeV and P — GeV^. 



the Minkowski space parton distributions through an inverse Melhn transform.^ More realistically, we might envisage 
extracting 6-10 even moments (and the parameters a and (3 in Eq. by fits using Eq. H18|l . Previous analysis |^ 
suggests that this is enough to reliably constrain standard parameterisations of parton distributions. Additionally, 
the low moments that have previously been computed directly from local operators can be used as input into our 
approach. In order to extract these higher moments, we must be careful of higher-twist contributions, as discussed 
in Sec. ^ Since by using a heavy quark we have removed many higher-twist terms, we might expect that those that 
remain will be small. This would be indicated by the fitted value of (3 being small. 

For simplicity, we work to zeroth order in the QCD coupling, taking the Wilson coefficients to be unit y. In a 
complete analysis, the perturbative Wilson coefficients appropriate to the desired scheme and scale [sH iBQllioL l4lj | 
should be used. With this assumption, the only unknowns in Eq. (|18(l are the moments ^^(/i) and the parameters 
a and /3. Then if we wish to determine six moments (n = 2, 4, . . . , 12) for example, we need the Compton amplitude 
evaluated at eight or more different combinations of momenta and heavy quark masses. A priori, we would also want 
to choose the masses and momenta such that the convergence of the expansion on the RHS of Eq. H18|l allows us to 
neglect terms beyond a certain n^ax and extract the moments for n < n^ax- For example, choosing = 3.54 GeV, 
go = 2.76 GeV and — 1.5 GeV^, the expansion on the RHS of Eq. ((TH|l falls off as shown in the left panel of Fig.^ 
However, from experimental measurements and perturbative counting rules, we know that the moments of the parton 
distributions in fact fall off rapidly as n increases and it may be more useful to choose the masses and momenta so 
that f{n) decreases slowly, allowing the natural suppression of the moments to control how many can be extracted. 



^ This is guaranteed to be a unique reconstruction by Carlson's theorem 1491 . Essentially we would be expanding the OPE in the Euclidean 
region, analytically continuing the Wilson coefficients and then re-summing the expansion in the Minkowski region. 
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FIG. 5: Extraction of moments of meson distribution amplitudes. Here, the light-shaded circle denotes the pion interpolating 
operator and the dark circle and dark square indicate the vector and axial- vector currents, respectively. 

In this case, choosing Af* = 2.1 GeV, qo = 1.98 GeV and = —3.85 GeV^ which gives the flatter behaviour shown 
in the RHS of Fig. ^ may be more appropriate. Without performing the large scale simulations that are required to 
determine the Compton amplitude, it is hard to be definite on the choices of parameters. However, it seems that this 
approach has significant potential to determine higher moments of isovector parton distributions than are currently 
available from QCD. 

IV. DISTRIBUTION AMPLITUDES FROM CURRENT-CURRENT MATRIX ELEMENTS 

A further application of the approach we have outlined is in computing moments of meson distribution amplitudes, 
In the lattice approach, we can extract moments of distribution amplitudes in the same way as DIS determines 
moments of parton distributions; for example, we may study the matrix element {tt^\T[V^ ^{x)A'^ ^{0)]\0) , where 
^ and A'^ ^ are fictitious vector and axial vector heavy-light currents. This process is described by the tensor 

= / d''xe^^-^7T+{p)\T[Vl^{x)Al^^,m\0) ■ (23) 

Following from Eq.@, the OPE of the two currents leads to the same matrix elements of twist- two operators that 
determine the moments of the pion distribution amplitude: 

{TT+ipM^^^^j^iiDr^ . . . {iOr'^HlO) = U{C~'). [P^' ■ ■ -P^" - traces] , (24) 

where 

(D.^ ['d^CMO- (25) 

Jo 

These matrix elements can be determined by studying the various components of S^"^^ for varying and 
analogously to Eq. 114|l . As in the DIS case, many higher-twist contributions are absent because of the valence nature 
heavy quark and the problems that plague direct evaluation of higher moments due to the lattice cutoff are eliminated. 
Since only the zeroth (deca y co nstant) and second moments of the pion distribution amplitude have been investigated 
in the direct approach [sil l52l IssL Isil Is^ . l56|| . any information on higher moments will be useful in constraining the 
distribution amplitude from QCD. For flavour non-diagonal mesons {e.g. tt^, K^'^), extraction of the tensor Sl^'^^^ on 
the lattice only requires the computation of the Wick contraction shown in Fig. 

V. SUMMARY 

To summarise, the direct study of Compton scattering tensor on the lattice using the operator product expansion 
can provide useful information on the moments of quark distributions. Using currents that couple an unphysical, 
quenched, heavy quark field to the physical light quarks renders the approach feasible without modifying the non- 
perturbative physics that can be extracted. This has the potential that a large enough number of moments can 
be extracted that the parton distributions can be reliably reconstructed from lattice calculations. Our analysis has 
focused on the unpolarised isovector quark distribution, but it can also be used to study the other twist-two and twist- 
three parton distributions and generalised parton distributions. Additionally, this method will allow computations of 
the moments of meson distribution amplitudes where even the lowest non-trivial moment is not known reliably from 
the lattice. 
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